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Abstract 

A phenomenological lattice model of high T c cuprates including order parame- 
ter phase fluctuations is considered within the BCS approximation, to interpret the 
experimental data from ARPES measurements on Bi2212 samples. A Kosterlitz- 
Thouless (KT) transition temperature T^ T is estimated below the mean field 
transition T^ F ', phase boundaries between competing order parameters of differ- 
ent symmetries are obtained and best model parameters, fitting the ARPES gap of 
d x 2_ y 2 symmetry, are determined. Variation of T^ T , as a function of the dopant 
concentration 5, is in qualitative agreement with experiments. 
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The symmetry of the superconducting order parameter (OP) in the high T c cuprate 
materials has been a widely debated issue in the last few years. In the recent past, 
several experiments done on the copper oxide materials [1-7], including those of phase 
sensitive experiments ||, interference measurements || and c-axis Josephson tunneling 
[[§], gave contrasting or inconclusive results regarding the OP symmetry. While a num- 
ber of experiments [p], |2j, || noted signatures of d x 2_ y 2 OP symmetry, others argued in 
favour of the anisotropic S"-wave or more exotic S + id wave OP symmetry [^, f|, |6], 
However, recently there has been considerable progress in the angle resolved photoemis- 
sion spectroscopy (ARPES) measurements || |9|], and a consensus seems to be emerging 
about the OP symmetry pi in the high T c cuprates. 



ARPES can give quantitative estimate of the momentum dependence of the super- 
conducting gap on the Fermi surface (FS) in terms of the spectral function representation 
of data and a detailed study of the FS is also possible [R [J. Due to its high angular as 
well as energy resolution || || , it can provide detailed and reliable knowledge about the 
nodes of the gap on the FS. Shen and coworkers, by ARPES measurements on Bi2212 
compounds M, found nodes of the superconducting gap on the FS along the 45 ° (it, tt) 
direction, suggestive of a d-wave OP symmetry. Ding et. a!., by similar measurements 
U on high quality Bi2212 single crystal, showed that the gap on the FS vanishes at two 
points (per quadrant) symmetrically displaced about the 45° direction, consistent with 
an anisotropic S-wave (S xy ) symmetry. However, reanalysis of their data revealed that 
the two node gap was an artifact of the superstructure producing "ghost" bands [[11 



In this light, the ARPES measurements on several B42212 samples were redone by Ding 
et. al. using dense sampling of the Brillouin zone (BZ) in the vicinity of the FS ||. The 
results are now consistent with a \cos(k x ) — cos(k y )\ type gap function implying a d x 2_ y 2 
OP symmetry. 

In order to interpret the two node gap data , a phenomenological BCS like lattice 
model in two dimension was introduced |T2| and mean field (MF) analysis of the instabil- 
ities in the spin singlet Cooper channel was done. The model is of interacting electrons 
on a square lattice, with an on-site repulsive (Vo) and attractive nearest neighbour (Vi) 
as well as next nearest neighbour (V2) interactions. The model Hamiltonian H = Hq + Hx 
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is 



and 



Hi = VoY^hijnu + Vy XX^+<W + V2 J2 J2 n ij)7 h i+ 5 nnn ,a> (lb) 



j, — * 

Where, el (cg CT ) is the quasiparticle creation (annihilation) opeator of momentum k and 
spin a, h is the number operator and 5 nn , 5 nnn are nearest neighbour and next nearest 
neighbour lattice vectors. The band dispersion eg was obtained by a six parameter 
tight binding fit to the normal state ARPES data on Bi2212 single crystal |12| where 



the parameters are [t ,...,t B ] = [0.131,-0.149,0.041,-0.013,-0.014,0.013] (in eV). Here 
to is the orbital energy, t\ nearest neighbour (nn), ti next nearest neighbour (nnn) etc. 
hopping matrix elements. Quasiparticle dispersion eg incorporates flat dispersion around 
(0, tt) and (n, 0) points which results van Hove singularity (vHS) in the single particle 
density of states (DOS). 



The mean field analysis of the model [0 found strong instabilities for the d x 2_ y 2 and 
S xy OPs which can best exploit the large single particle DOS just below the FS. Ratio 
of the interaction parameters V1/V2 determines the relative stability of these two states. 
Further extension of the mean field work, to include strong fluctuations present in the 
quasi two dimensional cuprate materials, was carried out and results were reported in a 
previous communication |13| . 



In this paper, we reanalyse the phenomenological model, including order parameter 
phase fluctuations, in view of the conclusive finding of a d x 2_ y 2 OP symmetry by the 
ARPES measurements on B42212 [|5], M. We first calculate a Kosterlitz-Thouless (KT) 
transition temperature T^ T , using the helicity modulus or superfluid phase stiffness 
(p s ) expression of the present model together with the KT relation p s {T^ T ) = ^kBT^ T , 
within each of the irreducible representation Bi (d x 2_ y 2), B 2 (d xy ) and Ai (S, S* and S xy ). 
We then find out the phase diagrams showing regions of relative stability of Bi, B2 and 
Ai states in the interaction parameter planes and determine the best model parameters 
corresponding to the d x 2_ y 2 gap with T^ T ~ 100 K. We also study the momentum 
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dependence of the gap function on the FS and dopant concentration (5) dependence of 
T^ T . Our main results are summarized below: 

(a) Best model parameters at the optimal doping level 5 = 0.17, with a Tf T ~ 100 K 
and for the stable B t (d x 2_ y 2) state, are V > 400meV, V\ ~ —48meV and 
V 2 85 meV. 

(b) Zero temperature d x 2_ y 2 gap on the FS, with V% = —48meV, match well the 
ARPES data |§ in the vicinity of 45° direction, but deviates from it as one moves 
beyond about 7 degrees on either sides. 

(c) With varying dopant concentration (8) , the gap magnitude at 0° angular direction 
changes considerably, although the node position is same for all 8 by virtue of the 
momentum dependent part of the d x 2_ y 2 gap function. 

(d) Tf T for B 1 (d x 2_ y 2) state function of 8 shows correct qualitative behaviour as 
in the high T c cuprate materials. 

Within the standard BCS approximation, the Hamiltonian of Eq.(l) yields the gap 
equation 

A £ = i E V(k - k>) ^ tanh (^j (2) 

where the quasiparticle energy is = J (eg — e^) 2 + | Ag| 2 and A^ is the BCS gap 
function. The pairing interaction, V(q) = V + 4Vi(cosq x + cosq y ) + 8V2 cosg x cosg y , 
is separable as V(k — k') = Ylt=o ViVi(k)Vi(k') ■ An expansion of the order parameter 
— Y^i=oVi(.k)Ai gives the linearized gap equation 

— — * 

^ = 4^^ tanh (^) E (3) 

where T) (k) = 1, 771 (A;) = ^(cosk x + cosk y ), r)2(k) = -(cos/c^ — cosk y ), T] 3 (k) = 
cos k x cos k y , 774 (fc) = sin k x sin k y are the basis functions corresponding to [S, S*, d x 2_ y 2, 
S xy , d xy ] symmetries, and (V0...V4) = (Vo, 8V\, 8V1, 8V2, 8V2). We suppress writing 
the terms corresponding to triplet pairing and ignore them in our analysis. Here 
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1Z = (Ai, Bi, B 2 ) are different irreducible representations of the C± v group and the 
gap equation factorizes to independent A2 (Bi representation), A4 (B2 representations) 
and three coupled linear equations A , A 1; A 3 (Aj representation). Since A 3 is the pre- 
dominant component within the Ai representation, it is identified as a state of S xy OP 
symmetry. 

To calculate the helicity modulus p s (and hence T^ T thereafter), a transverse vec- 
tor potential with the gauge A y = is considered. This introduces an extra phase 
which the carriers acquire while moving between the lattice sites. Hence the hop- 
ping matrix elements in H (Eq.(la)) should be changed through Peierls substitution 
tij — > Uj exp[|^ J 5* A ■ dl\. We work here with the units h = c = 1, but explicitly write 
them whenever necessary. 

The electron current operator j x (Ri) consists of the usual paramagnetic and diamag- 
netic terms [|14|]. To linear order in A XJ j x (Ri) is obtained by differentiating H with 
respect to A x (Ri) 

dA x (Ri) 

In Eq.(4), the paramagnetic term does not involve A x and is the electron velocity op- 
erator. The diamagnetic term is linear in A x and stems from the Meissner screening of 
the condensate. Average value of the diamagnetic current density is obtained as 

k,ir x 

where the ( ) represents an average in the mean field superconducting state. In a 
London like relation j x (q) oc — p s A x (q), the diamagnetic contribution to the phase stiff- 
ness (pf a ) is proportional to { c \ a c ka) ~dJ^~> ^ ne mean electronic kinetic energy 
along the x-direction JT5|| . The average in the lattice model turns out to be (cjr^Cfco-) = 
tanh ( : -f iL ) , unlike the continuum case where pf a oc — a (|^) ° 



e k~ e F 



(/ is the Fermi function). 



Contribution of the paramagnetic part is evaluated using linear response theory. In 
the long wavelength limit, the paramagnetic current is found to be 

1 



f x ara (q) = - 

c 



lim lim K xx (q,u) 



AM (6) 
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where K xx (q,u) = -i J dtO{t) e tujt {\jl ara {q,t), j% ara (-q,Q)}). The correlation function 
in Eq.(6) is evaluated to be K xx (q-> 0,cu -> 0) = jjy^Ejg (fjg) 2 ^M^- Taking the 
contributions from diamagnetic and paramagnetic parts, from Eqs.(5) and (6), we obtain 
the expression for superfluid phase stiffness 



k 



^V^ + ifell-^tanh^ 



dkj de % + 2dki y E n ' \ 2 



k 



(7) 

It should me mentioned here that, we work in a transverse gauge and vertex corrections 



required to get a gauge invariant current flIJJ have not been included. 



Above expression for p s involves linearized BCS gaps. In the inset of Fig.l we plot 
the superconducting gaps for different order parameters corresponding to B 1: B 2 and 
A\ states, as a function of temperature. The point where gaps become nonzero mark 
the mean field transition T^ F for a state. These gaps, as shown in the inset of Fig.l, 
are used to calculate p s {T) for different states. The KT transition temperature for 
each state is found by comparing the p s (T) curve for each state, with the KT relation 
p s (T^ T ) = ^T^ T . In Fig.l the intersecting point of a p s curve with the KT straight line 
(emerging from the origin) gives T^ T < T C MF . Similar technique was applied previously 
by Danteneer and coworkers for the two dimensional attractive Hubbard model 
and correct behaviour of T^ T , including its inverse coupling dependence in the strong 
coupling limit, was found. We too find a similar strong coupling dependence of Tf T 
in the present model. Thus, order parameter phase fluctuation degrades the mean field 
transition temperature. However, T^ T in our case is an upper bound of the actual KT 
transition temperature, that could be calculated only by considering superconducting 
gap renormalization due to the presence of vortex like fluctuations. 

Next, we consider the phase boundaries, calculated by comparing the T^ T values of 
the competing states. In Fig. 2, we plot phase boundaries separating d x 2_ y 2 and A\ states 
in the (-Vi, — Vz) plane, for various values of Vq. It clearly shows an widening of the the 
d x 2_ y 2 stable region with increasing on site repulsion for small values of —V\. Rate of 
this widening is faster for small Vq. For a fixed — V\, d x 2_ y 2 solution is preferred over A\ 
upto a maximum Vi- As for example, with Vq = 100 meV and V\ = —48meV, d x 2_ y 2 
solution is stable for — V2 < 80meV jT7]]. A comparison of the KT phase diagrams with 
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those from mean field calculations, can be found in Ref.[13]. 

To determine the optimal model parameters, one must also consider the phase bound- 
ary between A\ and d xy states and find out the region where A\ solution is not stable. 
In Fig. 3, we plot such phase boundaries in the (— V2, Vo) plane for different values of 
V\. To get the KT transition temperature T^ T ~ 100 K for d x 2_ y 2 OP, one must have 
V\ = —48meV. In Fig. 3, if we set V\ = —48meV and Vj> = —85meV, then to rule out 
the stability of Ax state, one must have Vq > 400 meV. Thus, we fix the optimal model 
parameters to be V\ = —A8meV, V2 = —85meV and Vq > 400 meV . 

With these optimal parameters, we study the momentum dependence of the zero 
temperature d x 2_ y 2 gap on the FS at 5 = 0.17. A plot of |A^(T = 0)| = |A 2 (T = 
0) x i[cos(fc 2; ) — cos(ky)} I is presented in Fig. 4 for two different values of V\. Solid circles 
are ARPES data || of T c ~ 87 K Bi2212 single crystal sample in the Y quadrant [0. We 
find that, our V\ = —48 meV curve matches well the experimental data till about 7° away 
from node position 45° on either sides, and falls below ARPES data beyond this. A much 
better match could be obtained with V\ = — 58meV. But, this makes T^ T ~ 135 K 
which is well above the sample T c ~ 87 K. One probable argument in favour of this could 
be that, our T^ T is an upper bound of the KT transition temperature. In actual case, 
gap renormalization due to the vortex fluctuations might reduce T£r . If one assumes a 
25% reduction, a new set of best model parameters can be obtained as V\ = —58 meV, 
V2 = — 90meV and Vq > 400 meV. In the inset of Fig.4, we plot the S dependence of 
the zero temperature gap which provides a feel for the d x 2_ y 2 gap magnitude at the FS 
angle = degree. 

Variation of T C XT for d x 2_ y 2 with dopant concentration 5 is presented in Fig. 5. This 
shows a correct qualitative bell shaped behaviour as in the high T c experiments. As 
a comparison, we also include Tf T for A\ state for the optimal parameters, which is 
below the d x 2_ y 2 T^ T curve for all 5. The T^ T curve here (and also A 2 vs 5 curve in the 
inset of Fig.4) peaks at around 6 = 0.25, that is consistent with the peak in the pairing 
density of states around e v ns ~ ~ 20meV which corresponds to 5 ~ 0.25. 

To summarize, we studied a phenomenological BCS model of high T c cuprate super- 
conductors including fluctuations, calculated phase boundaries separating order parame- 
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ters of different symmetries within the model and determined optimal model parameters 
to fit the d x 2_ y 2 superconducting gap as in the experiment. We also studied the dopant 
concentration dependence of the transition temperature which is in qualitative agree- 
ment with the experiments. 

Its a pleasure to acknowledge Dattu Gaitonde for useful discussions and suggesta- 
tions. 



8 



References 



[2 



[3 
[4 

[5 
[6 

[7: 

[8 
[9 

[10 

[11 
[12 

[13 
[14 



D. H. Wu et. ai, Phys. Rev. Lett. 70, 85 (1993); W. N. Hardy et. ai, Phys. Rev. 
Lett. 70, 3999 (1993). 

D. A. Wollman et. ai, Phys. Rev. Lett 71, 2134 (1993); ibid 74, 797 (1995); D. A. 
Brawner and H. R. Ott, Phys. Rev. B 50, 6530 (1994); C. C. Tsuei et. al, Phys. 
Rev. Lett. 73, 593 (1994); J. R. Kirtley et. ai, Nature (London) 373, 225 (1995). 

P. Chaudhari and S. -Y. Lin, Phys. Rev. Lett 72, 1084 (1994). 

A. G. Sun et. al. Phys. Rev. Lett. 72, 2267 (1994); R. Kleiner et. al. Phys. Rev. 
Lett. 76, 2161 (1996). 

Z. -X. Shen et. al, Phys. Rev. Lett. 70, 1553 (1993). 
H. Ding et. al, Phys. Rev. Lett. 74, 2784 (1995). 

S. Chakravarty, A. Sudb0 and P. W. Anderson, Science 261, 3999 (1993). 
M. Randeria et. al. Phys. Rev. Lett. 74, 4951 (1995). 



H. Ding et. al, Phys. Rev. Lett. 76, 1533 (1996); Preprint, |cond-mat/96"03044 
(1996). 

D. J. Scalapino, Physics Reports 250, 392 (1995); M. Sigrist and K. Ueda, Rev. 
Mod. Phys. 63, 239 (1991); J. Annet et. ai, Preprint, |cond-mat/960106q (1996). 

M. R. Norman et. al, Phys. Rev. B 52, 15107 (1995). 

R. Fehrenbacher and M. R. Norman, Phys. Rev. Lett. 74, 3884 (1995); M. R. 
Norman et. al, Phys. Rev. B 52, 615 (1995). 

B. Chattopadhyay, D. Gaitonde and A. Taraphder, Europhys. Lett, (in press). 

J. R. Schrieffer, Theory of Superconductivity, Addison- Wesley Publishing Com- 
pany, New York (1964). 



[15] D. J. Scalapino, S. R. White and S. C. Zhang, Phys. Rev. Lett. 68, 2830 (1992). 



9 



[16] P. J. H. Denteneer et. al., Europhys. Lett. 16, 5 (1991). 

[17] Below some values of —V\ and — V 2 , KT transition temperature T C XT < IK and 
its not possible to determine the phase boundaries beyond these parameter values. 

[18] In Y quadrant, the main CuO FS is widely separated from the two ghost FS and 
gap data are more reliable. Hence, Y quadrant data from Ref.[9] is chosen for 
comparison. 



10 



FIGURE CAPTIONS 



Fig.l. Helicity modulus p s is plotted as a function of temperature for different irreducible 
representations, at optimal doping 5 = 0.17 and with best model parameters Vq = 
400 me V, V± = —A8meV and V2 = — 85meV. The solid straight line, originating from 
(0,0) point, is the p s (T^ T ) = \T^ T line. [Inset: Mean field order parameters of different 
symmetries corresponding to A 1 , B\ and B 2 representations.] 

Fig. 2. Phase boundaries, separating the regions of stability of the A 1 and d x 2_ y 2 states, 
in the (— Vi, — V 2 ) plane, for 5 = 0.17 and various V values as shown in the figure. [Inset: 
Phase diagram for Vo = 0.4 eV is plotted on larger scales.] 

Fig. 3. Phase boundaries indicating the regions of stability of the A x and d xy states in the 
(— V 2 , V ) plane at 5 = 0.17 for various V\ values shown. Inset shows the phase diagram 
for V\ = —48 meV on an expanded scale. 

Fig. 4. Momentum dependence of the gap function |A k (T = 0)| on the Fermi surface is 
plotted against the Fermi surface angle 0, for 5 = 0.17 and for V 1 given in the figure. 
The angle is measured with respect to the line joining (n, n) and (0, n) points. Solid 
circles are ARPES data from Ref.[9]. [Inset: Zero temperature d x i_ y i gap as a function 
of 5.} 

Fig. 5. Kosterlitz-Thouless transition temperature Tf T for d x i_ y i and A 1 states as a 
function of dopant concentration 5, for optimal model parameters. 
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